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Abstract 

The cross section for prompt neutral and charged three pion production in pp inter- 
actions was measured at excess energies in the range 160 - 217 MeV. That comprises 
the first measurement of the pp — > ppTT°ir Tr° reaction and the direct comparison 
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with the pp — ► ppir + Tr~ir° process. The experiment was performed above the 77 me- 
son production threshold and the cross section could be directly normalized to the 
cross section of the pp — ► pprj reaction, with the 77 decaying into 3 pions. Since the 
same final states are selected, the measurement has a low systematical error. The 
measured cross section ratio a{pp — > ppTT + TT^TT°)/a(pp — > pptt°tt°tt ) is compared 
to predictions of dominance of different isobars in the intermediate state. 

ifey words: threshold measurement, three pion production, final state interaction 
PACS: 13.25.-k, 14.40.Aq, 25.40.Ve,25.75.Dw 



1 Introduction 



Production of three pions in proton-proton interactions, where the pions do 
not originate from decays of narrow meson resonances (like i] or a;), has 
not received proper attention yet, neither experimentally nor theoretically. 
Experimentally there are only very few data on the pp — > ppix+ix~ix° and 
pp — > pnit + -K + it~ reactions and no data on the other reaction channels: pp — > pp7T n n , 
pp — > pn-K + -K°it° and pp — > nn7i + 7Y + n . On the theory side there exists a com- 
plete microscopic model covering all reaction channels of two pion production 
in pN interactions [1], but so far no such models have been developed for the 
three pion case. In the isobar model the process pp — > pprnin should proceed 
by an excitation of one or two baryon resonances followed by the subsequent 
decays [2] and in the low energy region a mechanism with a simultaneous 
excitation of N* and A(1232)P33 resonances is expected to dominate. The 
iV* involved has to decay into Nirn and therefore the lowest lying Roper 
(JV(1440)Pn) and N(1520)D 13 resonances could be considered. There is how- 
ever a significant difference between their decay pattern. The Roper decays 
predominantly into An whereas for the iV(1520) the Np channel is equally 
important [3,4]. 

The influence of the resonances in the intermediate state can be studied in 
the invariant mass distributions of the subsystems of the outgoing protons and 
pions. Such studies were done for the pp — > ppn + 7i^7i and pp — > pn7i + 7i + n^ 
reactions in bubble chamber experiments performed at higher energies (beam 
kinetic energies of 4.15 GeV and 9.11 GeV) with up to thousand events [5-7]. 
In these studies, only one or two of the pions were considered to originate from 
decays of N* or A resonances. However, such kind of analysis is complicated 
and may not be conclusive due to many possible scenarios and due to large 
widths of the involved resonances. A simple starting point for analyzing the 
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three pion production close to threshold is to assume a constant value for the 
matrix element. The dependence of the cross section on the beam energy is 
then given by the phase space volume divided by a flux factor. For a proper 
description of the reactions in the threshold region it is required to take into 
account the final state interaction (FSI) between the outgoing protons [8,9]. 
The Ntt FSI on the other hand is expected to be negligibly small [10]. A hint 
about the production mechanism of the three pions could then be obtained 
by studying the ratio of the cross sections for the different charge states. The 
three pion production amplitude NN — > NNnnn can be written in terms 
of the isospin amplitudes M T . T37rT/ , where T { (Tf) denotes the initial (final) 
isospin of the nucleon pair and T 3n denotes the isospin of the produced pion 
triplet. It is now straightforward to show that 

aipp -> ppn+n-n ) oc ^|M 121 | 2 + ±|M m | 2 + ||M 101 | 2 

a( PP ^ PP 7l 7l 7l ) (x±\M in \ 2 

a(pp -> pnn+iT+ir-) oc ^|M 121 | 2 + ^|M m | 2 + ^|M 110 | 2 , 

where cross terms between Myi\ and M W i amplitudes are neglected. In the 
simple statistical approach [11], with all amplitudes M^r^Tf put equal, one 
obtains a{pp — > ppn + n^7Y°) : cr(pp — > ppn o n n ) : a{pp — > pn7Y + n + n^) = 8 : 
1 : 10. There were so far experimental data on a(pp — > pp7T + 7T~7T°) : a(pp — > pn7i + 7i + 7i^) 
only. At 2.0 GeV this ratio is 1:2.53±0.46 [12] and at 2.85 GeV it is 1:1.59±0.27 
[13]. This suggests a deviation from the statistical approach, especially at lower 
energies. Similar conclusions can also be drawn from experiments on double 
pion production at excess energies below 100 MeV [14] where the statistical 
approach fails. We discuss the cross section ratios under different assumptions 
about the dominating reaction mechanism in section 3 of this paper. 

With a 477 facility such as WASA [15], aiming for measurements of decays 
of i] and rj' mesons produced in pp interactions, the understanding of the 
pp — > ppmnr reactions becomes very important as they constitute a severe 
background for the studies of rj and rj' decays into three pions. Those decays 
provide key ingredients for determining the ratios of the light quark masses 
m u /m s and m^/m s [16,17], since the decay widths are proportional to the u 
and d quark mass difference squared. In order to estimate the effect of the back- 
ground for the decay experiments, the properties of the pp — > ppmiTf reactions 
for beam proton energies from 1.254 GeV (corresponding to the r] production 
threshold) to 3 GeV are important. There are no experimental points for the 
cross section of the pp — > pp-K°it -K° reaction in that energy range, and only 
three points for pp — > pp7i + 7i^7i° all coming from old bubble chamber experi- 
ments [12,13,18]. The experiment performed at the lowest energy, 1.48 GeV, 
identified a single event [18] (corresponding to 20 fib). Moreover in the bubble 
chamber experiments, the prompt three pion production was not separated 
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from the r] — > txtxtx decays. 



2 Measurement 

The analysis is based on data collected with the WASA facility at CELSIUS 
[19]. The target system provides small (0 ~ 30 \xm) hydrogen pellets that 
interact with the circulating proton beam. The protons have nominal kinetic 
energies of 1.30, 1.36 and 1.45 GeV (corresponding to excess energies, Q, in 
the center of mass system around 200 MeV for the pp — > ppirnir reactions). 
The integrated luminosities at each energy were: 27 nb _1 , 414 nb _1 (data on 
pp — > pp7T + 7T~7T° reaction are based on a sample corresponding to 80 nb^ 1 ) and 
221 nb^ 1 respectively. The WASA detector system consists of a multilayer 
forward detector (FD), for measurement of the outgoing protons scattered 
in an angular range of 2.5-18°, and a central detector (CD) containing an 
electromagnetic calorimeter and a drift chamber/solenoid for measuring the 
produced mesons and their decay particles in an angular range of 20-140°. 
The experimental method for extracting the cross section of the pp — ► pp-Kim 
reaction, relies on normalizing to the simultaneously measured pp — > pprj 
reaction, with subsequent decay of the eta into it + it~it (branching ratio 22.7% 
[4]) or into 7r°7r°7r° (32.5%). In this way, the reference reaction has the same 
particles in the final state, and most of the efficiency corrections cancel. The 
final state is selected by the requirement that two protons are detected in 
the forward detector. In addition the selected events should have six 7 hit 
clusters (with energy depositions of at least 20 MeV) in the calorimeter, for 
the 3"7r° case, or two 7 hit clusters and two tracks with opposite bending in 
the central detector drift chamber, for the 7r + 7r~7r° case. After identification 
of the two proton tracks in the FD no additional cuts on the proton-proton 
system are applied. The specific cuts for neutral and charged CD particles 
systems aim to select a clean 3tc final state. In case of the pp-K ir°iT channel, 
already the requirement of the six neutral hit clusters in the CD results in a 
fairly clean data sample. The 7r + 7r~7r° sample was obtained by requiring the 
invariant mass of the two photons to be located in the n° mass region, and 
the missing mass of the two protons plus the two photons to be greater than 
two pion masses. 

Two pion production, pp — > ppirir, is the major physical background which 
has to be considered in the analysis in view of its much higher cross section. 
The final state can be misidentified as a three pion event only in case of two 
additional hit clusters in the calorimeter, due to cluster split-offs, pile-up or 
noise. To reject the background strict time cuts are applied, as well as addi- 
tional kinematical conditions - such as reconstruction of the total energy. The 
final contribution of the background is on the per cent level, in agreement with 
expectations from the Monte Carlo simulations. Therefore the uncertainty in 
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Fig. 1. Experimental pp missing mass distributions for ppmnr final states with fit 
of the Monte Carlo distributions for 2tt, 3tt and r\ production. 

the exact value of the total cross section for the two pion channels, as well as 
in the reaction mechanism, only gives a minor contribution to the systematical 
error. 



The derivation of the cross sections for the pp — > pp3ir reaction is based on 
the fact that the resolution in the pp missing mass determination is very good 
(4 MeV/c 2 FWHM at 1.30 GeV) and that the missing mass distribution is 
practically insensitive to the mechanism of the pp — > ppr] reaction. 

After the selection of the ppimn final states, the pp missing mass distribution is 
constructed from experimental data, and from Monte Carlo simulated data of 
the reaction channels pp — > pp-innr, pp — > pp(i] — > 3ir), and pp — > ppnir as the 
main background contribution. The individual Monte Carlo distributions were 
then fitted to the measured spectrum, with the pp — > ppnnn to pp — > ppr] ratio 
as parameter. The pp — > ppnn contribution was fixed relative to the pp — > pprj 
channel based on the cross section values listed in table 1. Fig. 1 shows the 
result of the best fit at the three beam energies for the pp — > pp7i 7i°7i° reaction 
and at 1.36 GeV for the pp — > pp7i + 7i^n reaction. Several studies were done 
in order to determine uncertainties and check consistency: 

• For 1.36 GeV, data from two well separated run periods were analyzed and 
gave consistent results. 

• The kinematical cuts were varied. 

• Application of a kinematical fit to the pp — > ppmrii reaction and verification 
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Fig. 2. Comparison of experimental distributions in the lab system of a) proton 
kinetic energies, b) 7 reaction angles, and c) pp opening angle for beam proton 
kinetic energy 1.36 GeV projectiles with Monte Carlo simulation of the pp3ir° , pp2ir° 
and pp(rj — > 3ir°) best fit including pp FSI. 

of the x 2 distribution (7 constrains (7C) fit for the tt tt tt° and 5C fit for 
the 7T + 7r~7r°). 

• Variation of the nominal beam energy and of the energy resolution of the 
detectors assumed in the analysis. 

The applied method to obtain the cross section ratio relies heavily on the 
Monte Carlo simulations. Their accuracy is demonstrated in fig. 2 by compar- 
ing distributions of some variables with the result of the simulation, taking 
into account all three mentioned reaction channels. The normalization factors 
for the three channels are the same for each distribution, and were obtained 
from the fit of the pp missing mass distribution (fig. 1). The simulation de- 
scribes well detector effects such as the structure at 0.3 GeV, which is due to 
an ambiguity in distinguishing stopped protons from punch through protons 
in the FD. The angular distribution of the 7s shows two distinct, well repro- 
duced structures. Near 40° the binning changes to account for different sizes 
of the Csl scintillators in the CD; at 90° the detection efficiency is reduced by 
the pellet target components. The major geometric acceptance limitation is 
imposed by the angular acceptance for protons (3°-17°). The good agreement 
between simulated and real data both in the gamma and proton scattering 
angle distribution indicates that the acceptance correction is under good con- 



6 




.70 



M 2 pp [(GeV/c 2 ) 2 ] 



Fig. 3. Proton-proton invariant mass squared at 1.36 GeV. The solid line shows the 
parametrization from this experiment. The dashed line shows the model prediction 
from A. Deloff [20] tuned to describe the TOF data [21] and dotted line the pp FSI. 

trol. 

In order to extract the total cross section, a model has to be used to extrap- 
olate data outside of the measured range of the proton scattering angles. For 
the pp — > ppr] reaction the largest deviation from the phase space behavior is 
expected due to the strong pp final state interaction. We use a full calculation 
of the pp FSI of A. Deloff [20] with the Reid A/W-potential, which yields a 
perfect reproduction of the 1 Sq phase shifts for relative pp momenta up to 
0.3 GeV/c. This calculation reproduces the pp angular distribution much bet- 
ter, but not yet the squared missing mass (MjL) distributions obtained for 
Q = 15.5 MeV [22] and Q = 41 MeV [21]. Agreement was then obtained 
in [20] by expanding the production amplitude in the r\ momentum relative to 
the pp pair. The proportionality constant was optimized for a best fit of the 
Mp p distribution, cf. fig. 3. Application of this approach to the M pp distribu- 
tion of this work yields a very similar result for Q — 41 MeV (see fig. 3) and 
Q — 17 MeV (not shown). This fit at the same time gives a good description 
of the pp opening angle distribution (fig. 2c), where the pp FSI leads to an 
enhancement at the lowest angles. For the pp — > pp3n, the outgoing protons 
can have scattering angles up to 30° at those beam energies. This means that 
extraction of the total cross section involves some extrapolation into unmea- 
sured kinematical regions. However, it was checked that e.g. the inclusion of 
a broad resonance in the Ntt system in the reaction mechanism changes the 
total acceptance by a few per cent only. 



3 Results and discussion 

The main result of the experiment is the ratio between the pp — > pptttttt and 
pp — > ppr] cross sections for which values are given in table 2 with statistical 
and systematical uncertainties. One can extract the total cross section values 
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Table 1 

Total cross sections for two pion and eta production used in the analysis. The data 
for 7r°7r° production are extracted from bubble chamber experiments [18, 23] and 
from CELSIUS/WASA [24]. The data for vr+vr" production are from ref. [25]. The 
data for eta production are from CELSIUS, COSY and Saclay [22,26-28]. Some 
values were obtained by interpolation to the energies used in this experiment. 



Beam kinetic 


PP 7 ? ppvr + 7r 


pp7r°7r° 


energy [GeV] 


[fib] [fib] 


[fib] 


1.30 


2.64 ±0.25 


(1.6±0.4)xl0 2 


1.36 


4.9±1.1 660±100 


(2.0±0.3)xl0 2 


1.45 


16±2 


(3.5±0.8)xl0 2 



Table 2 

The results for the ratio of the cross sections of the pp — > ppTnnr reaction and the 
pp — > ppr] reaction with the corresponding r\ decay into 37r. Both statistical and 
systematical uncertainties are shown. In addition the extracted cross sections for 
pp — > pp7r TT°TT° and pp — > ppir + ir~Tr° reactions are listed. 



Beam kinetic 


pp7r 7r°7r° 


a(pp - 


■+ ppTT°TT TT°) 


energy [MeV] 


ratio prompt/77 




[fib] 


1299±2 


0.480±0.053±0.11 


0.42 ±0.05 ±0.10 


1361±2 


0.550±0.014±g;g3 


0.89 ±0.02 ±0.23 


1448±3 


0.260±0.010l8;Si 




l±0.05t°j° 


Beam kinetic 


ppir + ir~ir° 


a{pp —> 


ppir + 7r~ir°) 


energy [MeV] 


ratio prompt/77 




[fib] 


1361±2 


4.1±0.3±0.4 


4.6 ±0.3 ±1.2 



using the known a ppv from table 1 together with the relevant 77 decay branching 
ratio into 3ir. For the beam proton energy at 1.30 GeV the 77 cross section 
points come from four experiments [22,26-28]. At 1.36 GeV there is only one 
data point [26] to compare with. The cross section value at 1.45 GeV, is an 
interpolation using the data from PINOT [27]. The extracted values of cross 
sections for the production of the three pions are also given in table 2. 

As a cross check of the acceptance corrections for the pp — > ppr] reaction, the 
cross sections have also been estimated using the luminosity derived from the 
simultaneously measured pp elastic scattering events [29]. The scattering cross 
sections were taken from the precision experiment EDDA [30]. The obtained 
values for the three energies 1.30, 1.36 and 1.45 GeV are: 3.36±0.17±0.5 /xb, 
5.06±0.26±0.6 fib and 14.9±0.75±1.0 fib respectively. The uncertainties in- 
clude statistics, luminosity determination (5%) and the acceptance correction 
(including uncertainty of the reaction model). We find a general agreement 
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Fig. 4. Measurements of the pp — > ppTr + TT^ir° and pp — > ppTr°Tr ir° cross section in 
the beam proton energy range 1—3 GeV. The four experimental points below 1.5 
GeV are from the present experiment, while the data at 2 GeV and at 2.85 GeV for 
the pp — > ppir+Tr~ir° reaction are from bubble chamber experiments [12, 13]. Dotted 
lines are predictions of the statistical model [11] (phase space and flux factor). 
Dashed and solid line predictions are with pp FSI from Deloff [20] and from Faldt 
and Wilkin [8] respectively. 



with the a PPV results from literature, which indicates that detection accep- 
tances are well under control. 

Figure 4 summarizes the results for the reactions pp — > pp7i + n^n° and pp — > ppn n n . 
The solid lines include the pp FSI calculated according to the parameteriza- 
tion from Faldt and Wilkin [8] and the dashed lines the pp FSI calculated 
with the Reid wave function [20] (here, for pp relative momenta in center of 
mass greater than 300 MeV/c, pure phase space is used). The lines from the 
model predictions are normalized to the experimental data points at 1.36 GeV. 
In addition the dotted lines give the energy dependence of the cross section 
calculated from the statistical model and are normalized at high energies to 
the dashed lines. Two data points for pp — > pp7Y + n^n° from bubble chamber 
experiments in the energy range 2-3 GeV [12, 13] are also shown. The data 
point at 2 GeV [12] has been corrected by subtracting the fraction of events 
expected from the pp — > ppr] reaction [31]. 

The experimental value obtained for the ratio of the cross sections for pp — > ppii + 7i^7i 
and pp — > ppir -K°-K° at 1.36 GeV is 5.2±0.5(stat)±0.8(syst). To have a compar- 
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ison at the same Q, the result should be corrected for the difference between 
the charged and neutral pion masses, which makes the phase space volume at 
1.36 GeV for the pp — > ppn°n°n° reaction 18% larger than for pp — > ppit + it~it c '. 
The threshold corrected value for the ratio a(pp — > ppir + n~n ) /a(pp — > ppir°n°n°) 
thus becomes 6.3±0.6±1.0. As mentioned in the introduction a similar ratio 
is expected in the statistical approach with all amplitudes being equal. The 
simplest two baryon excitation mechanism involves the formation of an inter- 
mediate A and an N* (1440), where the iV*(1440) decays either to N (7r7r)Jr° ave 
or via an intermediate An state [4]. In the iV*(1440) — > N(nn)J=® ave scenario 
T 37r has to be I. Then only the amplitude M m can contribute and a cross 
section ratio of 4 is expected from the expressions given in section 1. This 
decay branch can thus certainly be a major part of the total reaction mecha- 
nism. For the scenario iV*(1440) — > An, T = 0,1 or 2 is permitted, and thus 
all M 101 , Mm and M m amplitudes can be involved and the ratio for the 
3n system can not be calculated without further assumptions. Note however 
that in this case the experimental ratio can easily be reproduced by choos- 
ing certain values of the amplitudes Mioi and Mm- For example the choice 
|M 101 | 2 Rs 0.7|M m | 2 and |M m | 2 = gives the observed ratio 6.3. Accordingly, 
also the reaction diagram involving iV*(1440) — > An might well be the leading 
part of the reaction mechanism. This conclusion is further supported by cal- 
culations within the isobar model by Sternheimer and Lindenbaum [2] which 
lead to the ratio 7. In case of 7V(1520) in the intermediate state one expects 
a much larger value, since the additional N p decay mode can contribute only 
to the pp — > ppn + n~n° channel. 



Comparison of the presented result with the cross sections obtained from 
the Crystal Ball data [32] on n~p — > 37r°n allows to rule out a mecha- 
nism with sequential iV(1535) decays. In both pp and n~p interactions the 
i] production is dominated by the excitation of the iV(1535) resonance and 
close to the thresholds one expects that o(n~p — > 37r°n) jo{n~p — > rjn) ~ 
a{pp — * pp?>n°) / a{pp — > pprj) for this scenario. However, the ratio a{ii~p — > 
3n°n)/o-(n~p — > rjn) was measured to be below one percent [32], much lower 
than the 10% for a{pp — > pp3n°)/cr(pp — > pprj) presented in this paper. 



In conclusion microscopic model calculations, of the same kind as those exist- 
ing for the double pion production [1], are needed to shed more light on the 
issue and exploit the result of the presented measurements. Experimentally, 
more data are desired to get information on the cross sections for the other 
three pion reaction channels. 
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